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Abstract 

In proposals for quantum computers using arrays of trapped ultracold polar molecules as qubits, a strong 
external field with appreciable gradient is imposed in order to prevent quenching of the dipole moments 
by rotation and to distinguish among the qubit sites. That field induces the molecular dipoles to undergo 
pendular oscillations, which markedly affect the qubit states and the dipole-dipole interaction. We evaluate 
entanglement of the pendular qubit states for two linear dipoles, characterized by pairwise concurrence, as 
a function of the molecular dipole moment and rotational constant, strengths of the external field and the 
dipole-dipole coupling, and ambient temperature. We also evaluate a key frequency shift. Aw, produced by 
the dipole-dipole interaction. Under conditions envisioned for the proposed quantum computers, both the 
concurrence and Atj become very small for the ground eigenstate. In principle, such weak entanglement 
can be sufficient for operation of logic gates, provided the resolution is high enough to detect the Acli shift 
unambiguously. In practice, however, for many candidate polar molecules it appears a challenging task to 
attain adequate resolution. Simple approximate formulas fitted to our numerical results are provided from 
which the concurrence and Aw shift can be obtained in terms of unitless reduced variables. 
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I. INTRODUCTION 



Since the original proposal by DeMille UJ], arrays of ultracold (< 1 mK) polar molecules have 
come to be considered among the most promising platforms to implement a quantum computer [O- 
13n . His proposal describes a complete scheme for quantum computing using as qubits the dipole 



moments of diatomic molecules, trapped in a one-dimensional optical lattice, partially oriented 
in an extemal electric field, and coupled by the dipole-dipole interaction. The qubit states are 
individually addressable because the field has an appreciable gradient so the Stark effect is different 
for each location in the array. 

A subsequent proposal has advocated coupling polar molecules into a quantum circuit using 
superconducting wires [|l4ll . Such capacitive, electrodynamic coupling to transmission line res- 



1611 . The molecular 



onators is analogous to coupling to Rydberg atoms and Cooper pair boxes [|15L 
qubits are entangled via the coupling to the transmission lines rather than direct dipole-dipole in- 
teractions. Again, addressability of the qubits is achieved via the Stark effect by means of local 
gating of an electrostatic field. 



Entanglement is a major ingredient in most quantum computation 



defining features of quantum mechanics, with no classical analog [Il7l-ll9h. A pure state of a pair 



algorithms. It is among the 



of quantum systems is said to be entangled if its wavefunction cannot be factored into a product 
of wavefunctions of the individual partners. For example, the singlet state of two spin-^ particles, 
■^(1 Ti) - I iT)) is entangled. A mixed state is entangled if it cannot be represented as a mixture 
of factorizable pure states. The allure of quantum information processin g ha s recently motivated 



20I- I32I1 . These include 



studies of entanglement for a variety of potential qubit systems p, la, |9|, 
one-dimensional arrays of localized spins, coupled through exchange interactions and subject to 
an extemal magnetic field [l23ll and analogous treatments of trapped electric dipoles coupled by 



dipole-dipole interactions [|32ll . 



However, the previous studies of entanglement of electric dipoles have not adequately consid- 
ered how the extemal electric field, integral to current designs for quantum computers using polar 
molecules, affects both the qubit states and the dipole-dipole interaction. For the simplest case of a 

diatomic molecule, the qubit eigenstates resulting from the Stark effect are linear combinations 
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of spherical harmonics, with coefficients that depend markedly on the field strength. These are ap 



propriately termed pendular states [13311 . or field-dressed states [|34ll . In such states, the orientation 
of the dipole moment has a broad angular range (not solely along or opposed to the field direction 
as are spins in a magnetic field). Likewise, the dipole-dipole interaction for molecules in pendular 
states is much different than that for dipoles in the absence of an external field. 

Here we evaluate entanglement, as measured by pairwise concurrence, for the prototype case 
of two diatomic polar molecules in pendular states, ultracold and trapped in distinct optical lattice 
sites. The molecules are represented as identical rigid dipoles, undergoing angular oscillations, 
a fixed distance apart and subject either to a different or to the same external electric field. We 
examine the dependence of the concurrence on three dimensionless variables. The first governs 
the energy and intrinsic angular shape of the qubits (when the dipole-dipole interaction is switched 
off). It is /i£/5, the ratio of the Stark energy (magnitude of permanent dipole moment times 
electric field strength) to the rotational constant (proportional to inverse of molecular moment of 
inertia). The second variable governs the magnitude of the dipole-dipole coupling. It is Q/5, with 
Q = (jU^/r^), the square of the permanent dipole moment divided by the cube of the separation 
distance. The third variable, ksT/B, is the ratio of thermal energy (Boltzmann constant times 
Kelvin temperature) to the rotational constant. 

We also examine an aspect related to but distinct from entanglement. The operation of a quan- 
tum gate I35I1 such as CNOT requires that manipulation of one qubit (target) depends on the state 
of another qubit (control). This is characterized by the shift, Aco, in the frequency for transition 
between the target qubit states when the control qubit state is changed. The shift Aco, which is due 
to the dipole-dipole interaction, must be kept smaller than the differences required to distinguish 



m for 



among addresses of qubit sites. Under conditions envisaged in the proposed designs [|1|- 
quantum computing with trapped polar molecules, Q./B < 10 and for the ground eigenstate both 
the entanglement and frequency shift Aoj become very small. For CNOT and other operations, en- 
tanglement needs to be large, but can be induced dynamically, so need not be appreciable in the 
ground eigenstate. Yet a small Aa* shift can only suffice if the resolution is high enough to detect 
the shift unambiguously. From estimates of the line widths of transitions between the pendular 
qubit states, we find it an open question whether adequate resolution can be obtained for typical 



candidate diatomic molecules. 



II. ENTANGLEMENT FOR TWO DIPOLES IN PENDULAR STATES 
A. Hamiltonian terms and pendular qubit states 

The Hamiltonian for a single trapped linear polar molecule in an extemal electric field is 

2 

Yi = ^ + Vrrap{r) + Bf-n-Z (1) 
2m 

where the molecule, with mass m, rotational constant B and body-fixed dipole moment //, has 
translational kinetic energy /Im, potential energy Vnup within the trapping field, and rotational 
energy as well as interaction energy /i-£ with the extemal field £. In the trapping well, at 
ultracold temperatures, the translational motion of the molecule is quite modest and very nearly 
harmonic; I2m + Vtrap{v) thus is nearly constant and can be omitted from the Hamiltonian. There 
remains the rotational kinetic energy and Stark interaction, 

Hs = Bf - pEcose (2) 

which represent a spherical pendulum with 6 the polar angle between the molecular axis and the 



field direction. Figure HI a) displays the lowest few pendular eigenenergies [13 611 for a diatomic 
(or linear) molecule, as functions of pE/B. These are labeled with the familiar quantum numbers 
/, M that specify the field-free rotational states. However, / wears a tilde to indicate it is no 
longer a good quantum number since the Stark interaction mixes the rotational states, whereas M 
(denoting the projection of the J-vector on the field direction) remains good as long as azimuthal 
symmetry about £ is maintained. As proposed by DeMille, the qubit states |0) and |1) are chosen 
as the lowest M = pendular states, with / = and 1, respectively. These are superpositions of 
Yjfl spherical harmonics, 

10) = ^jYjM 1 1 ) = 2 bjYjM (3) 

j j 

Figure |2] plots the coefficients as functions of pEjB. Figure |3] displays the angular distributions 
of the pendular qubit states. For |0) the distribution is unimodal and as pEjB increases the dipole 



orientation increasingly favors the direction of the £-field (at 6 = 0"). For |1) the distribution 
is bimodal because, with M = 0, the dipole is rotating perpendicular to the J-vector, which is 
perpendicular to the field direction. For £ = 0, the dipole orientation is equally probable in the 
hemispheres toward (6 < 90") or opposite (6 > 90") to the field direction. As ju£/B increases, the 
pinwheeling dipole favors the opposite hemisphere because there its motion is slowed because the 
Stark interaction becomes unfavorable. However, when luE/B becomes large enough, pinwheeling 
is inhibited and converted into pendular libration about the field direction, so the dipole orientation 
shifts to favor the toward hemisphere. 

Adding a second trapped polar molecule, identical to the first but distance apart, introduces 
in addition to its pendular term the dipole-dipole coupling interaction, 

w /fi H2-3(pi •n)(//2 -n) 
|ri - r2p 

Here n denotes a unit vector along ri2. In the presence of an external field, it becomes appropriate 
to express Vd-d in terms of angles related to the field direction. As shown in Appendix A, the result 
after averaging over azimuthal angles (that for M = states are uniformly distributed) reduces to 

Vd-d = ^(1 - 3COS^Q')COS01COS02 (5) 

where Q.= jjL^ /r^^, the angle a is between the Yn vector and the field direction and polar angles 6\ 
and 6*2 are between the IX\ and H2 dipoles and the field direction. Until later (Sec. IV), we consider 
the extemal field magnitude and direction to be the same at the sites of both the polar molecules. 



B. Entanglement measured by pairwise concurrence 



We will deal with the entanglement of formation, (£(p), which characterizes the amount of 
entanglement needed in order to prepare a state described by a density matrix, p. (Henceforth, we 



term (£(p) just "entanglement", for short.) Wootters [|37 



3811 has shown that C(p) for a general state 



of two qubits can be quantified by the pairwise concurrence, C(p), which ranges between zero and 
unity. The relation can be written as 

(£(p) = ^(C(p)) (6) 
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where f is given by 

^(C) = h 



(7) 



2 

with h{x) = -x\og2X - (1 - x)log2(l - x). The function ^(C) increases monotonically between zero 
and unity as C varies from to unity. The concurrence is given by 

C{p) = max (O, - - - V^j (8) 

where the /i,'s are the eigenvalues, in decreasing order, of the non-Hermitian matrix pp, where p is 
the density matrix of the spin-flipped state, defined as 

P = {ay®(Ty)p*{(Ty®(Ty) (9) 

with p* the complex conjugate of p and o-y a Pauli matrix. The parent density matrix p is taken 
in the basis formed by combining the pendular qubit states; for a pair of two-level particles, this 
comprises the four state vectors {|00), |01), |10), |11)}. 

In order to evaluate thermal entanglement, we need a temperature dependent density matrix, 
p = exp(-/3H)/Z(T), with /3 = l/Z^^r and Z(T) the partition function 

Z{T) = tr[cxp{-m] = gie-"""- (10) 

i 

with Ei the eigenvalue and gi its degeneracy. Hence the density matrix can be written as 

1 ^ 

P(r) = -2e-^^'|¥,)<^,| (11) 



where I*?,) is the z"' eigenfunction. From the density matrix p{T), we can obtain the reduced 
density matrix for any pair of dipoles and thence evaluate the concurrence at any temperature. 

III. CONCURRENCE OF TWO DIPOLES IN PENDULAR STATES 

We illustrate the calculation of pairwise concurrence for N = 2 dipoles. The Hamiltonian, Hsi 
+ Hs2 + Vd~d, when set up in a basis of the qubit pendular states, {|00), |01), |10), |1 1)}, takes the 
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form, 



Hsi + H 



S2 



Wo + VF^ 

Wo + W[ 
1^1 + 












Wi + W[ 



(12) 



Vd-d = ^^(1 - 3cos%) 



(13) 



CqC'q CqC^ ^xC'q CxC'^ 
CqC'^ CqC\ CxC'^ CxC[ 
CxC'q CxC'^ ^1^0 ^i^x 
CxC'^ CxC[ CiC^ CiC[ 
where Wo and Wi are the eigenenergies of the pendular qubit states |0), and in the absence of 
the dipole-dipole interaction. Primes attached to quantities for the second dipole indicate that the 
external field magnitude may difl'er at its site (although, as noted above, we postpone evaluating 
that case until Sec. IV). In Vd~d the basis qubit states are linked by matrix elements containing 
factors arising from the orientation cosines in Eq. (|5}; these are 



Co = <O|cos0|O); Cx = <O|cos0|l); Ci = <l|cos0|l) 



(14) 



Co and Ci are the expectation values of cos6' in the pendular states |0) and |1), respectively, so 
represent for those states the eff'ective dipole moment projections displayed in Fig{TJb). Cx cor- 
responds to an exchange interaction or transition dipole moment between the qubit states. Both 
the Stark eigenenergies Wi and the dipole-dipole elements Ck are functions of /uE/B. As seen in 
FigHtb), as fxE/B is increased Co becomes increasingly positive, whereas Ci is increasingly neg- 
ative until about /J.2/B = 2, then climbs to zero at about fi£/B = 4.9 and thereafter is increasingly 
ositive. The range jjE/B = 2 to 4 is recommended for the proposed quantum computer designs 



po; 

& 



1611 : within that range, the diff"erence in the effective dipole moments of the qubits, |Co - Ci|, 



varies only modestly. 

If the dipole-dipole interaction is omitted (Q. = 0), the eigenvectors of Hsi + Hsi are simply 
^1 = 100), = 2-^'^(\l0) - 101)), = 2-1/2(110) + 101)), = |11), corresponding to the 
eigenenergies of Eq.(12). For ^Fi and ¥4, which are obviously nonentangled states, the concur- 
rence is zero. For *I'2 and ¥3, which exemplify fully entangled states, the concurrence is unity; 



these are termed Bell states nlln . 

When the dipole-dipole coupling is included, an analytical solution to obtain eigenstates is only 
feasible when the external field is switched ofl". As shown in Appendix B, in that limit analytical 
results can be obtained for each step in evaluating the concurrence, both for the four individual 
eigenstates and their combination in the thermal concurrence. As seen in Fig. 1, for fi£/B = 0, the 
energy terms in Eq.(12) involve merely Wq = and Wi = IB. In the Vd-d matrix of Eq.(13), the 
cosine matrix elements Co and Ci vanish and Cx = 3"^''^; thus, the only nonzero elements occur 
along the antidiagonal and (for a = 90") are just QC|. The results for this zero-field limit prove 
useful in interpreting those for the general pendular case. 

The limits with ^1 = and/or /i£/5 = motivate setting up the Hamiltonian of Eqs.(12) and 
(13), for the (unprimed) case with the same external field at both dipole sites, using a basis of Bell 
states: 

111) + 100) 111) -100) 110) + 101) 110) -101) ^^^^ 



V2 V2 
In this basis, the Hamiltonian becomes 
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■2 



Where = W^+Wo and = ^(C^+CD+Cl, B = ^(C^-C^), 4 = Cx(Ci+Co), = CiCq+C 
This makes explicit a consequence of the symmetry between the (unprimed) sites [|39ll . In the Bell 
basis, the Hamiltonian factors, with the state I'^^^Q 10) - 101 )) in a 1 x 1 block, so that state remains 
maximally entangled regardless of the value of fi£/B or Q/B. 

Figure IDplots, for fiE/B = 0, 2 and 4.9, the eigenenergy and pairwise concurrence versus Q./B 
= to 6 for the four eigenstates of the two-dipole system. The eigenstates are numbered from 
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1 to 4 in order of increasing energy. For /i£/5 = 0, both eigenstates 2 and 3 are Bell states, 
with eigenenergies EJE = 2 - (Q./6B) and 2 + {Q./6B), respectively; eigenstates 1 and 4 are 
also entangled (much more weakly) by the dipole-dipole interaction, with eigenenenergies that 
shift downwards and upwards nonlinearly with increasing Q./B, respectively. For //£/5 > 0, the 
concurrences increase with Q./B for eigenstates 1 and 4, and decrease for eigenstate 3. By virtue of 
the symmetry imposed factorization noted above, eigenstate 2 retains the same Bell form despite 
the Stark and dipole-dipole interactions which affect its energy, and its concurrence is always unity. 
For small Q./B « 1, eigenstate 3 also becomes independent of the dipole-dipole interaction 
and coincides with eigenstate 2 in both energy and concurrence. For fi£/B = 4.9, as seen in 
Fig. 1(b), the Ci = (l|cos6'|l) factor that appears in seven of the matrix elements in Eq. (13) 
vanishes. Consequently, the energy of eigenstate 4 then becomes independent of the dipole-dipole 
interaction, although its wavefunction and concurrence do not. 

Figure |5] shows, for fiE/B = and 2, how the contributions of the basis states to each of the 
eigenstates vary with the strength of the dipole-dipole interaction. This illustrates that for Q./B « 
1 the eigenstates rapidly approach those for Q = 0. Indeed, we find that for Q/B < 0.04 the 
concurrences for eigenstates 1 and 4, which rapidly become the same, are proportional to Q./B 
within better than 1%. Thus, 

Ci2 = K(x)[a/B] (18) 

where the proportionality factor K{x) is a function of jc = fi£/B. At the zero-field limit, ^(0) = 1 /6. 
In Appendix C we describe a numerical analysis that provided an accurate approximate formula, 

K(x) = Ai + — (19) 

1 -I- exp[(x - xo)/Ax] 

This is plotted in Fig. 6 and values of the four parameters are listed in Appendix C. 

Figure |7]displays for Q./B = 0.1, 1 and 6 the eigenenergies and concurrences versus fi£/B from 
to 8 for the four eigenstates. As the dipole-dipole interaction increases 60-fold over this range, its 
effect on the eigenstate energies is relatively modest, whereas the concurrences change markedly, 
in response to variations in eigenvector compositions such as illustrated in Fig. |5l 

Figure [8] gives a contour plot of the thermal pairwise concurrence derived from Eq.(ll) as a 
function of Q./B and ksT/B. It pertains to juE/B = 3; we found that normalizing the thermal 



concurrence to its value for T = and D./B = \ removed most of the variation with /i£/5 from 
such contour plots. For T = 0, the thermal concurrence coincides with that for the ground state, 
eigenstate *Fi. However, as ksT/B increases, the thermal concurrence decreases and is always 
smaller than the ground-state concurrence. This may seem odd, because Eq.(l 1) specifies a shift in 
population that reduces the contribution from the gound state, while bringing in contributions from 
the excited states. The eigenstates 2 and 3 then populated have large concurrence, so increasing 
temperature might be expected to make the net thermal concurrence become larger than for the 
ground-state, rather than smaller. The source of this behavior is indicated by the analytic solution 
obtained in Appendix B for the zero-field limit, 

Ci2(T) = Cn(l)Pi - Cu(2)P2 - Cn(3)P3 - Cu(4)P, (20) 

where Pi = {\ /Z)exp(-Ei /kBT) with Z(T) = Jjiexp{-Ei/kBT). This shows that the excited states 
indeed reduce the thermal concurrence, an effect traceable to Eq.(8) and which persists even for 
large /uE/B. 

Another striking aspect of Fig. 8 is that the concurrence vanishes along and outside a particular 
contour. That contour defines mutually dependent maximum values of ksT/B and minimum values 
of Q./B required to obtain nonzero concurrence. When Q./B « 1, we find that a modified form of 
Eq.(fT8l) represents the thermal concurrence, 

Ci2(T) = max{0, Kix)[y - yo{x, z)]} (21) 

Here x = fi£/B; y = Q./B; and z = k^T /B is the scaled temperature. Fig. |9]gives a contour plot 
of Jo = ^c/B, the critical dipole-dipole coupling required for nonzero concurrence. Some further 
details are included in Appendix B. 
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by DeMille and kindred papers on quantum computing with trapped polar 
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16D discuss for several examples the range of experimental conditions 
deemed suitable and acceptable. For trap temperatures of the order of a microkelvin or below, 
the typical values of ksT/B are a few times 10"^, so indicate that only ground-state entanglement 
would be significant. The extemal field strengths considered are typically a few kV/cm. The 
spacing between optical lattice sites, r = A/2, is half the optical lattice wavelength. The optimal 

10 



choice of A ranges between 1 to 0.3 microns, depending on electronic transition frequencies of 
the molecules to be trapped ll^. From these parameters and molecular data, values of Q./B are 
small; we find for a dozen potential candidate molecules values ranging between 4 x 10"^ (for 
KCs) to 2 X lO""^ (for Csl). A favorite candidate is SrO (jj. = 8.9D, B = 0.33 cm-\ A = 1 micron), 
for which Q./B ~ 10^. In that regime, the concurrence is simply proportional to Q./B, so can be 
easily evaluated from Eq. (fT8] ) and/or (|2T|) without use of the rather elaborate prescription outlined 
in Eqs.(6-14). 



IV. FREQUENCY SHIFT FOR TWO COUPLED DIPOLES IN PENDULAR STATES 

In the region Q/B < 10 the concurrence of the ground eigenstate is very small, typically 
< 10"^. However, such meager entanglement in eigenstates can still be adequate for quantum 



computing, as demonstrated with NMR versions of quantum computers [|40|] . The key aspect is that 
although entanglement needs to be large for some quantum computing algorithms, it need not be 
appreciable or even present in the ground eigenstate of the system; it can be induced dynamically 
during operation of the computer [l4lil . Here, for the polar molecule case, we consider this aspect. 
We also evaluate an eigenstate property, a small frequency shift, distinct from but related to the 
pairwise concurrence, that is important for quantum computing. 



430 is an essential 



The need for selective excitation in operation of quantum logic gates [|35L 
feature. Taking the 2-qubit CNOT gate as an example, its operation requires that manipulation of 
one qubit (target) is perceptively affected by the state of the other qubit (control). In our case, the 
qubits are pendular states that can be accessed by microwave transitions, which offer high spectral 
resolution. As resolution has a crucial role, we now suppose the external field differs enough at 
the two dipole sites (denoted unprimed and primed) to supply distinct addresses for the sites (cf. 
Fig. 1(a), green dashed curve). 

Since Q./B is so small, we first omit the dipole-dipole interaction and, as illustrated in Fig. [TOl 
consider transitions among the pendular eigenstates of Eq.(12). Although in this limit the ground- 
state concurrence is zero, as seen in Eq. (fr8T ). it is possible to generate states of large concurrence 



by use of resonant pulses [|4lL 



44|] . Start by applying a pulse resonant with the transition denoted 
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ioy, between |00) and |01), which has energy W[ - Wq. Note that (jJ\ needs to be well-resolved from 
the transition ojj,, between |00) and |10), which has energy Wy - Wq. The separation thus comes 
from the different values of the external field at the two sites (plus a dipole-dipole contribution, 
in higher order). The requisite field strength difference, £' - £, can be readily determined from 
another approximation formula, 

iWy - Wo)/B = Ai + f] (22) 
1 + {x/xo)' 

by comparing results for x = luE/B and x' = ju'E'/B; the accurate fit obtained (better than 1% except 
near x = 0) is displayed in Fig. [TT]and the four parameters in Eq.(E21i are given in Appendix C. 
The amplitude and duration of the oji pulse can be adjusted to make it a n/2 pulse, which will put 
the system in the state 2"'^^(|00) + |01)). 

Next, to complete the CNOT gate, apply a pulse resonant with the transition 0)2 between |01) 
and 111). This needs to be well-resolved from transition to^ between |00) and |10). However, in 
our initial approximation, both t02 and 0*3 have the same transition energy, V^i - V^o- Hence, weak 
as it is, the dipole-dipole interaction is seen to have an essential role: to introduce a frequency 
shift, Aio = iOs - CO2, adequate for unambiguous resolution. If that is fulfilled, the amplitude and 
duration of the 0*2 pulse can be adjusted to make it a ;r pulse. Thereby the system will be put in 
the state 2"^^^(|00) -1- 11 1)). This result of a CNOT gate is to first approximation a Bell state (aside 
from small corrections of order Q./B), so its concurrence will be near unity. It is not an eigenstate, 
so will evolve with time but in principle would remain nearly fully entangled until degraded by 
other interactions. 

If now the dipole-dipole terms from Eq.(13) are included to first order, we obtain 
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(23d) 



where = - 3cos^Q'). Thus, the key frequency shift is given by 

t\(X) = (X)^- ix)2 = (x)^- (x)i= ^.^iCi - Co)(C'i - Cg) (24) 
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For given Q.^, the frequency shift Ato depends only on x and x', which determine at the respective 
sites the difference in the effective dipole moment projections Q and Ci along the extemal electric 
field, specified in Eq.(14). To provide a convenient means to evaluate Eqs.(23) and (24) we again 
fitted our numerical results to obtain accurate approximation formulas, 

Co{x)=A,+ (25) 
1 + {x/xoY 

Ci(jc) = Ao + — + — (26) 

1 + exp[(jc - Xi)/Axi] 1 + exp[-(x - X2)/AX2] 

These functions are plotted in Fig. [TTl together with Co - Ci, and the fitted parameters are given 
in Appendix C. 

Since for small Q./B, both the concurrence and Aco are proportional to Q./B, the frequency shift 
provides an equivalent measure of entanglement. When the £-fields differ at the two sites, Eq.(18) 
still provides a very accurate approximation for Cuix, x'), merely by replacing the proportionality 
factor by the geometric mean, [K(x)K{x')y^^. The concurrence (which involves Cx, the exchange 
interaction term) is in principle different from Ao; but both have about the same magnitude. The 
frequency shift is much more relevant for quantum computing, because Aco is directly involved in 
the CNOT gate. 

Also important, in addition to the pulse shapes which affect the population transfers, are the 
durations of the resonant pulses required to resolve a)\ and 0)2 from 0*3; these must satisfy T31 >> 
1/|cl»3 - coil and T32 >> 1/|<^3 - a»2|. For T31 the lower bound usually can be made very low, 
permitting a short pulse duration. This holds because a£ as well as dipole-dipole terms contribute 
to \coj - a»i|, which thus can be made large by choice of the £-field gradient, regardless of whether 
Q.a is extremely small. In contrast, for T32 the separation Aoj = |cl»3 - 0*2 1 depends only on the 
dipole-dipole interaction. The smaller Aco is, the longer the 0J2 pulse duration has to be in order 
to complete the CNOT operation. Although larger Aco allows a shorter pulse duration, Aco must 
not be so large that it becomes comparable to or larger than the addressing shift produced by a£, 
thereby thwarting correct identification of the qubits. 
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TABLE I: Illustrative CNOT Gate Implementation''. 





x= 1 


x' - 1.01 


x' = 1.10 


X - 3 


x' = 3.03 


x' = 3.30 


{Wi - Wo)/B 


2.2709 


2.2759 


2.3218 


3.5614 


3.5831 


3.7789 


Co 


0.30165 


0.30404 


0.32487 


0.57922 


0.58149 


0.60051 


Ci 


-0.16467 


-0.16573 


-0.17461 


-0.16362 


-0.16150 


-0.14115 


Co - Ci 


0.46632 


0.46977 


0.49948 


0.74284 


0.74298 


0.74165 


(a>i - W3)/S 




4.99E-3 


5.09E-2 




2.17E-2 


2.17E-1 


Aco/B 




2.19E-6 


2.33E-6 




5.52E-6 


5.51E-6 


Cn 




1.20E-6 


1.17E-6 




3.57E-7 


3.34E-7 



"For Qa/B - 10 ^; x - jj£,/B, x' - jj^' jB. As the quantities in the lowest three rows are functions of 
both X and x' , their values are listed in the x' columns. There E-n denotes a factor of 10 ". 



Table I provides specific numbers pertaining to the SrO example. From Sec.III, we take ^alB = 
10'^. As representative £-field values, we use x = 1 and 3 for site 1 and and take x' higher by 1% 
or 10% for site 2. From Eqs.(23), the transition frequencies coi = 0)4 and 0*2 = 0*3 (in units of B) 
to 5 or 6 significant figures. The frequency difference that must be resolvable for the first step of 
the CNOT operation, oji - coj,, is approximately just AAW = (W[ - W'^) - (Wi - Wq). From Fig. 
[TTl this is seen to grow about linearly with both x and x' - x. The values in Table I (third row from 
bottom) range from > 10"^ to > 10^' (in units of B). To accommodate more dipole qubits, it may 
be desired to make much smaller the £-field differences between sites; steps with Ax = 0.01% 
were proposed by DeMille [uJ]. That might encounter engineering limitations, but in principle 
the proportionally smaller cji - 0)3 difference could still be readily resolved. For the second step 
of the CONT operation the crucial frequency shift, Aoj = oji - toi, varies only modestly with x 
and practically not at all with x' - x. The values of Aco/B in Table I (second row from bottom) 
range between 2 and 6 x 10"^; in frequency units, this range is 20 to 60 kHz. Smaller still are 
the corresponding values of the concurrence (bottom row); also insensitive to x' - x but, in accord 
with Fig. 6, varying more rapidly with x. 

Figure [l2]exhibits for both the Aa> shift and concurrence the variation with a, the angle between 
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the direction of the electric field and the axis between the dipoles. This dependence enters via the 
factor (1 - Scos^or) in the dipole-dipole interaction, Eq.(5), which emerges directly in the Aa» shift, 
Eq.(24), and by a more complex route propagates into the concurrence, via Eq.(13). Tilting the 
field direction to make a = 54.73°, the "magic angle", provides a siniple means to shut off the 



entanglement. That is a useful option, awkward to attain in other ways [|7, 



m. 



V. CONCLUSIONS AND PROSPECTS 

In this study, our chief aim has been to examine entanglement of polar molecules by the dipole- 
dipole interaction and subject to an external electric field, for the prototype case of two diatomic 
or linear ^ 2 molecules. This required use of qubits that are pendular states comprised of sums of 
spherical harmonics. We focused on the pairwise concurrence and its dependence on three unitless 
reduced variables, involving the dipole moments, field strength, rotational constant, dipole-dipole 
coupling and temperature. We have considered a wide range of the parameters, to map general 
features of the concurrence. However, for conditions envisioned for proposed quantum computers, 
the dipole-dipole coupling is weak {Q./B typically of order 10""^ to 10"^) and the concurrence 
becomes very small (< 10^^). For that weak coupling realm, we found the Aoj frequency shift 
provides an equivalent measure of entanglement, directly related to observable properties and 
hence preferable to the concurrence. We also obtained for both the Acj shift and concurrence in 
the weak realm simple explicit formulas in terms of the reduced variables. 

For quantum computing a crucial issue is whether Aoj is large enough to enable the toj transition 
to be reliably distinguished from 0)3 (and, equivalently, coi from C04). For typical candidate polar 
molecules, this requires resolving transitions separated by only tens of kHz. That would not be 
feasible in conventional molecular spectroscopy. Under ordinary gas phase conditions, transitions 
between molecular rotational or pendular states have line widths of the order of a few 100 kHz 



[|45n . For ultracold molecules trapped in an optical lattice, line widths may be much narrower. 
CoUisional broadening is eliminated and at microkelvin temperatures Doppler broadening is also 
quenched (as trap conditions are in the Lamb-Dicke regime). It is encouraging that for ultracold 
atoms extremely narrow line widths have been attained by exploiting "magic" optical trapping 
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conditions that are expected to be at least in part applicable to molecules 14611 . At present, however, 
no data have been reported on line widths for rotational transitions of ultracold molecules trapped 
in an optical lattice and subject to an external electric field. In view of the small size of Aa», it 
is important to obtain such data to assess the resolution attainable, since motion within the traps, 
coupling to lattice fields, and inhomogeneity of the external field may introduce appreciable line 
broadening. 

We have sought to glean pertinent evidence from electric resonance spectroscopy of molecular 
beams, as the beams are collision free and transitions are observed in an external electric field 
("Rabi C-field"). For BaO, both AM = 0, 7 = ^ 1 transitions in the microwave region [Lvj and 



Ay = 0, \M\ = 0^1 transitions in the radiofrequency region [|48|| have been observed, in fields 
ranging from ~ 200 - 500 V/cm. For the radiofrequency transitions, line widths were only about 
2 kHz, consistent with just the dwell time in the C-field. But for the microwave transitions the 
widths are much larger, 45 kHz; this is attributed both to the higher frequency of the transitions 
and to experimental conditions that render more significant Doppler broadening and nonunifor- 



mity of the field, especially in the entrance and exit fringe regions [|49|1 . The Doppler and dwell 
time contributions are not relevant to inferring what might be expected for trapped BaO (or SrO). 
Broadening by inhomogeneity of the external field is relevant but depends very much on experi- 
mental particulars. The transitions of interest, depicted in Fig. \Wi occur in the microwave region 
and involve Stark fields typically ten-fold larger than used in the electric resonance spectroscopy, 
so the line widths might be significantly broadened due to field inhomogeneity. These observations 
do not permit firm conclusions about the resolution issue, but it decidedly poses an experimental 
challenge. 

This discussion pertains only to the choice of qubits we have considered, pendular states of 
linear polar molecules, which involve transitions that change / but not M. The resolution issue 
motivates examining other choices for qubits. For instance, states with the same / but different 
M could be used. Other options, particularly use of hyperfine or nuclear spin states instead of 
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m- As 



pendular states, have been suggested as means to reduce sources of decoherence yj, 
yet, the size of Aco for any qubit choice other than that used in this paper remains to be determined. 
We intend to extend the treatment developed here to other choices for qubit basis states as 
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well as to larger numbers N > 2 of dipoles. In preliminary work on linear and planar arrays of 
dipoles up to A'^ = 8, we find, as expected, the maximum pairwise concurrence occurs for next- 
neighbor dipoles, although that for non-nearest ones is significant. Also in prospect is an analogous 



treatment of the proposed coupling of polar molecules via microwave strip-lines u6n . There the 
entangling interaction differs from direct dipole-dipole interactions, but again the proposed qubits 
are pendular states. 

Another pendular variant inviting attention is use of polar symmetric top molecules. The |0) 
and |1) qubits can be selected as |7, ^, M) = |1, 1, -1) and |1, 1, -l-l), which are degenerate in the 



field-free limit and thus have a first-order Stark effect [14511 . Even in a weak electric field, these 
states are strongly oriented along and opposed to the field, with equal and opposite projections. 
Moreover, the effective dipole moments do not depend on the field strength so low fields can be 
used if necessary to reduce line broadening, without the penalty imposed by quenching of effective 
dipoles that would occur for the second-order Stark effect. At first blush, the symmetric top option 
appears to be disallowed because transitions between the M = -I and -l-l Stark states violate the 
selection rules, AM = or AM = +1. But the prohibition is not absolute. Because the optical 
lattice perturbs cylindrical symmetry about the field, M is not strictly a "good" quantum number, 
so the AM selection rule is relaxed. Moreover, if the molecule contains an atom with nuclear spin 
I > 1/2, and hence an electric quadrupole moment, transitions with AM = +2 become allowed. 



For instance, a deuterium nucleius (/ = 1) makes AM = ±2 transitions facile in Stark spectra [|50||. 
Other symmetric top options for qubits are inversion doublets (e.g. in NH3) or internal rotation 
states associated with hindered torsional motion (e.g. CH3CF3); these offer strong dipole- allowed 
transitions. 

Previous studies of entanglement, both for polar molecules [Q, 



[123 
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32|1 and for magnetic spins 



31[1 . have considered primarily domains where the concurrence is large (> 0.1), and have 
focused on means to tune the entanglement to attain such domains. For polar molecules, that 
requires Q./B > 1. Recently, it was suggested that such large Q./B could be attained for dipole 
arrays by exploiting nanotraps with lattice spacing of the order of only 10 nm [|32t1 . However, 
as emphasized in Sec. IV, for quantum computing large entanglement in the ground eigenstate 
is not required. Indeed, reducing the array spacing so markedly would strongly foster inelastic, 
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spontaneous Raman scattering of lattice photons and hence induce decoherence [Hi |2|, HSfl. Such 
considerations make small rather than large and consequently weak rather than strong en- 
tanglement in the ground eigenstate, actually preferable for quantum computing [|37ll . provided 
resolution of the ho) shift can be attained. 
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APPENDIX A: DIPOLE-DIPOLE INTERACTION 



The angular dependence of the dipole-dipole interaction, given in Eq.©, is usually expressed 

as 

0,-y = coSjS - 3cosy;Cosy^- (Al) 

where ji is the angle between dipoles /i, and //y; angles y, and yy specify the orientation of the 
dipoles with respect to the vector r,-, between them. Ordinarily, it is natural (and done in all 
textbooks) to express coSjS in terms of the angles y together with the azimuthal angles (pr about the 
r,y axis. Thus, use 

coSyS = cosy,cosyy -I- sinyisinyycos(0r; - <Prj) (A2) 



which when combined with the -3cosy,cosyy term gives the familiar expression [|51ll . In the pres- 
ence of the extemal electric field, we need to recast 0,y in terms of angles Oi and 6j that specify the 
orientation of the dipoles with respect to the direction of the extemal electric field. Therefore, we 
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use 

cos/3 = cosOiCOsOj + s'm6iS'm6jCOs{(f>Ei - 4>Ej) (A3) 

cosy = cos^/coscr + sin0;sinacos(0£,- - (f>Er) (A4) 

where the (pE are azimuthal angles about the field vector E and a is the angle between the field 
vector and the r,y vector. The azimuthal factors can be expressed as 

cos(0i7,- - (pEj) = cos(f)EiCOS(pEj + s'mcpEiSm^Ej (A5) 

As M = states, which do not depend on the (p angles, are chosen as the qubit basis states, in 
evaluating matrix elements of O/y between these states the integrations over dcpidcpj (from to 2n) 
eliminate all terms involving the (pE angles. The net result is simply 

(^ij)^ = (1 - 3cos^a)cos0,cos6'^- (A6) 

The effect of integrating over the 6 angles is just to replace in Vd^i the dipole moments and yUy by 
their effective values, ju{cos6). The effective dipole-dipole interaction hence reduces to 

Vd-d = ^^(1 - Scos^aXcos^/Xcos^y) (A7) 

with Q. = fXiiUj/r^j as a convenient scale factor. 



APPENDIX B: ZERO-FIELD CASE 

For yu£/5 = 0, the Hamiltonian matrix reduces to diagonal terms from Eq.(12) and antidiagonal 
elements from Eq.(13), and the pendular qubit basis states become simply |0) = Foo and |1) = Yw 
The form of the Hamiltonian makes it equivalent to that for the Ising model for a system with two 



qubits 1521] . Diagonalization of the Hamiltonian yields the eigenenergies and eigenvectors given 



in Table n as explicit functions ofQ/B. 
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TABLE II: Zero-Field limit for N = 2 dipoles. 



/ Eigenenergies, EjlB 



Wavefunction, 4', 



1 


2 - 2(1 + 


2 


2(1-0 


3 


2(1 + 


4 


2 + 2(1 + 



(|ll>-a+|00» 



^(I10>-|01» 
^(|10> + |01» 
=(|ll>-a_|00» 



/f, with (^^l/6B 



where a± - 1 + ( 1 + 
The density matrix for eigenstate 1, the ground-state, is 



p(l) = l^iX^il 



l+ai 








-a+ 




















-a+ 





1 



That for eigenstate 2 is 



^0 0^ 



p(2) = |»P2)(^2l = +2 



C 



12 



2a+ 

1 
1 

1+0-1 



(Bl) 



(B2) 



1-10 
0-110 
^0 0^ 

The p(3) matrix differs from p(2) by having -l-l in place of each -1; the p(4) matrix differs from 
p(l) by having a- in place of a+. As these density matrices pertain to only two dipoles, they need 
not be reduced further. 

Obtaining the density matrices, p(z), for the spin-flipped states, defined in Eq. (|9l), involves 
shuffling the rows and columns of p{i) in accord with |00) [11) and |01) |10). This gives 



P(l) = 



l+ai 



' \ -an ^ 






-Q'+ y 



(B3) 
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and thus the product matrix is 



P(1)P(1) = 



la\ -a\ 






-la^ lal 



(B4) 



+ / 



The eigenvalues of this matrix are X\ = Aa\l {l + a^j , = h = M = From Eq. ([8]), the 
concurrence is, Cnil) = 2a+/{l+al_). Similarly, we find the concurrences for the other eigenstates, 
given in Table I. 

To evaluate the thermal concurrence, we need to set up the thermal density matrix 

/ 



where 



with 



Then we find 



4 



p(T) = ^oxp(-Ei/ksT) l^'tX^t 



i=l 



a g 
b d 
d b 
[g c) 



l+ai 



1 +a^_ 



c = 



b = -(Pi + P3) 



\+ai 



1 +a^_ 



1 



d=-iP3-P2) 



0^4 



p(T)p{T) 



l+ai I + at 

_ expi-Ej/ksT) 
Z 

ac + g^ 
b^ + d^ 2bd 
Ibd b^ + d^ 
leg 



lag 



ac + g^ 



(B5) 



(B6) 

(B7) 
(B8) 

(B9) 
(BIO) 

(BID 



(B12) 
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and obtain the eigenvalues, 

Ai=(^^ac-gf, A2 = (b-df, A^ = {b + df, A^ = {^fac + gf (B13) 

Hence from Eq. ([8]), we obtain the themial concurrence, 

CniT) = max {0, -2{b + g)} (B 14) 

With the Cnii) of Table I, this gives Eq.(20) of the text, 

Cn{T) = CniDPi - Cn(2)P2 - Cn(3)P3 - Cn(4)P4 (B15) 

When Q./B « 1, the ground state concurrence becomes Ci2(l) Q./6B, in accord with 
Eq.CEl) of the text, whereas Ci2(4) ^ n/6B and Cuil) = CuiS) = 1. Provided that also 
kgT IB « 1, a first-order expansion of Eq.(B15) gives 

Ci2(r) « Ci2(l) - 62 - 63, (B16) 

where 62 = PilPi « 1; 63 = P3,IP2 « 1 and P^/Pi «< 1. Then 

CniT)^K{0)[n/B-QJB], (B17) 
which has the form of Eq.(2 1 ) of the text with K{0) = 1/ 6 and 

3;o(r) = QJB = (62 + €3)/K{0) = Uexp {-2B/kBT) cosh (Q/SkBT) (BIS) 

This result for the zero-field case, although not useful in practice, illustrates how the excited 
states are involved in creating a temperature dependent minimum level of dipole-dipole coupling, 
O-c/B, that is required to have nonzero thermal concurrence. 

Figure [12] shows a contour plot of C^iT) for the zero-field case, derived from Eq.(B14). It is 
qualitatively quite similar to Fig. 8 for the pendular case, over the same wide range of ksT/B and 

n/B. 

APPENDIX C: REDUCED VARIABLE FORMULAS 

In order to find a proper reduced variable formula, three steps are needed: (1) calculate enough 
sample points to define well the exact curve; (2) find a function with adjustable parameters that 
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enables fitting those points; (3) evaluate the parameters using a non-linear regression method. 



For our curve fitting we use the Levenberg-Marquardt Algorithm 115 3h . also called "Chi-square 
minimization". Chi-square is defined as: 

N 

x^ = J][Yi-fixr,e)f (CD 

!=1 

where Xj and F, are the independent and dependent variables for the i''^ (i = 1, 2,...,n) sample points 
of the exact curve; 6 are the parameters to be fitted. The Levenberg-Marquardt algorithm iteratively 
adjusts the parameters to get the minimum chi-square value, which corresponds to the best fit. The 
input data for fitting Eqs.(19, 22, 25, 26) comprised our numerical results for the pendular case, 
over the ranges x = juE/B = to 8. Tables III - VI list the optimal values found for the parameters 
and 95% percent confidence intervals. At x = 0, the Eq.(19) fit gives ^(0) = 0.17103, diff"erent 
slightly from the exact zero-field limit, K(0) = 1 /6. Likewise, at x = the Eq.(25) and (26) fits 
give Co = 0.005 and Ci = 0.00072 rather than the exact value of zero. The critical point for Ci to 
change sign is x = 4.902, whereas Eq.(26) gives Ci = -0.00025, slightly different from zero. 

TABLE III: Values of the parameters for Eq.([T9]l. 



Paiameters 


Pendulai- 


Field-free 


CF 


Ai 


0.01092 


0.00221 


+0.0003 


A2 


0.21953 


0.24779 


±0.006 


xo 


0.96578 


0.74035 


±0.05 


Ax 


0.97429 


0.86072 


±0.03 



"95% confidence interval; values listed are maximum found for the 2 curves 
shown in Fig|6] Both R~ values are around 0.9981. Similarly accurate results 
are found when Eq.lfTsTl is generalized for different E-fields at the dipole sites 
by replacing K{x) by [K{x)K{x')y'-. 
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TABLE IV: Values of the parameters for Eq.(l22]i. 



Parameters Values CF' 

Ai 12.42379 ±0.0533 

A2 -10.47646 +0.0534 

;co 8.77516 +0.0534 

p 1.5867 ±0.00527 



" 95% confidence interval. = 0.9999 



TABLE V: Values of the parameters for Eq.(l25Tl. 



Parameters Values CP 

Ai 0.84855 ±0.00145 

A2 -0.84355 ±0.00180 

;co 1.6339 ±0.00508 

p 1.2459 ±0.00539 



" 95% confidence interval. = 0.99994 



TABLE VL Values of the parameters for Eq.(|26]l. 



Parameters Values CI' 



Ao 


-0.75212 


±0.0323 


Ai 


1.04192 


±0.0336 


A2 


1.14092 


±0.0325 


XI 


-0.16241 


±0.0224 


X2 


3.1232 


±0.124 


Axi 


0.90544 


±0.0136 


AX2 


2.76286 


±0.0496 



" 95% confidence interval. R- - I 
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For convenience, we give formulas for the three unitless ratios, evaluated with customary units: 

//8/5 = O.Omiu(Debye)£(kV/cm)/B(cm~^y, 

Q/B = 5.04 X \Q-'^ti\Debye)lr\microns)IB{cm~^); 

hT/B = 0.695 T{K)IB{cm-^). 
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FIG. 1: (Color online) Stark states for a polar diatomic molecule in a electronic state Il36il . as functions 
of yu£/B, with jU the permanent dipole moment, £ the field strength, B the rotational constant, (a) Eigenen- 
ergies, W, and (b) Matrix elements of orientation cosines; see Eq.(14). States used as qubits (red curves) 
are labeled |0) and In the field-free limit, |0) correlates with the 7 = 0, My = and |1) with the 7 = 1, 
M/ = rotational states. Dashed curve (green) in (a) shows energy for transition between qubit states, 
= Wi - VKo; that in (b) shows Co - Ci, difference between eff'ective dipole moments, projections of the 
molecular dipole on the field direction for pendular states |0) and |I). 
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FIG. 2: (Color online) Coefficients of sums of spherical harmonics for pendular states |0) and |1>, see Eq.(3). 
Dashed curve for |1) indicates the coefficient of yo,o is negative. 



29 



30 60 90 120 150 180 30 60 90 120 150 180 

^ (degrees) ^ (degrees) 

lo 



FIG. 3: (Color online) Angular distributions of the |0) and |1) pendular states for values of juE/B between 
and 8. 
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FIG. 4: (Color online) Eigenenergies, numbered 1, 2, 3, 4 in order of increasing energy, and pairwise 
concurrences of the eigenstates for two dipoles as a function of the dipole-dipole coupling constant HjB for 
three values of the reduced electric field strength, fi£/B = 0, 2 and 4.9. 
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FIG. 5: (Color online) Eigenvectors of the four eigenstates for two dipoles as a function of the dipole-dipole 
coupling constant D./B = to 6 for /uE/B = (dashed curves) and 2 (solid curves). 
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FIG. 6: (Color online) The K{x) function of Eq.(18); solid curves show exact result (red) and fitted function 
(blue) of Eq.(19) that pertains to pendular qubit basis, see Eq.(3). For comparison, dashed curve pertains to 
field-free basis with |0) = 7o,o, |1> - J'l.o (Cf. Table III, Appendix C). 
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FIG. 7: (Color online) Eigenenergies and concurrences for the four eigenstates for two dipoles as a function 
of reduced variables, iu£/B for electric field and Q./B for dipole-dipole coupling. 
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FIG. 8: (Color online) Contour plot of thermal paiwise concuixence for two dipoles, for iu£/B = 3. For 
< Q./B < 1, the maximum concurrence Cnimax) = 0.0473, occurs at T = 0, Q./B = 1. The plot displays 
normalized contours. Within each colored band, the variation of Ci2/Ci2{max) is 0.1; thus the nomialized 
concuiTcnce in the right most band (red) ranges from 0.9 to 1 ; in the next band (orange red), from 0.8 to 0.9, 
etc. A striking feature is the large region (uncolored) where Cu = 0. There, entanglement does not occur 
unless the dipole-dipole couphng exceeds a critical value dependent on the temperature. 
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FIG. 9: (Color online) Contour plot displaying yQ{x,z) term in Eq.(21) vs. x = ij£/B. Within each colored 
band, the range of yo is 0.001; thus in the lowest colored band (magenta). 3^0 ranges between and 0.001; 
in the highest colored band (red), 3^0 is between 0.009 and 0.01. 
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FIG. 10: (Color online) Schematic energy levels for qubit pendular eigenstates of N = 2 dipoles, in absence 
of dipole-dipole interaction, thus corresponding to Eq.(12). Qubit basis states shown at left, eigenenergies 
at right. Pairs of transitions involved in CNOT operation are indicated: coi transfers dipole 2 from |0) to 
|1> with dipole 1 remaining in |0>; then a)2 transfers dipole 1 from |0) to |1) with dipole 2 remaining in 
|1). Analogously, the same result could be reached by ojt, followed by 0)4. Transition energies (including 
dipole-dipole terms to first-order) are given in Eqs.(23). 
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FIG. 11: (Color online) Comparison of exact results (blue curves) with fitted approximation functions 
(dashed red cui^ves) for properties governing transitions among qubit states, Eqs.(23): pendulai^ energy 
diff"erence, {Wi - Wo)/B, cosine expectation values, Co and Ci and their difference, Co - Ci; cf. Eqs. (|22l 

USES. 
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FIG. 12: (Color online) Frequency shift Aa>/B (left panels) and ground state concurrence C12 (right panels) 
as functions of a, the orientation angle of the electric field. Cui^ves are shown for Q./B = 0.2 to 1.0 with 
fiE/B = 1 or 3. 
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FIG. 13: (Color online) Contour plot of themial paiwise concurrence for field-free case, prepared in same 
format for comparison with Fig. 8 for the pendular case. Here, Cuimax) = 0.1644 at T = 0, Q./B = 1. 
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